function [gestureID,likelihood] = test_sequence(sequence,model)

    alpha = 0.001;

    sequence = {getMatrix(sequence)}; % Convert to a 2D matrix
    
    load('model.mat','model');
    
    for j=1:length(model)
        % Compute LL of kth sequence of class i using model j
        loglik(j) = mhmm_logprob(sequence,...
            model{j}.hmm.prior,...
            model{j}.hmm.transmat,...
            model{j}.hmm.mu,...
            model{j}.hmm.sigma,...
            model{j}.hmm.mixmat);
    end

    % Scale by the length of the sequence (frames)
    loglik = loglik/size(sequence,2); 

    [loglikMAX, gestureID] = max(loglik);
    
    % compute confidences of each model j for its
    % predicted label from (scaled) log likelihoods
    % P(x|y)
    Likelihood = exp( (loglik) .* alpha );

    % Confidence of predicted label
    likelihood = Likelihood(gestureID);

end